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Abstract 

Two particles that are just shy of binding may develop an infinite number of shallow bound 
states when a third particle is added. This counterintuitive effect was first predicted by V. Efimov 
for identical bosons interacting with a short-range pair-wise potential. The Efimov effect persists 
for non-identical particles, if at least two of the three bonds are almost bound. The Efimov effect 
has recently been verified experimentally using ultra-cold atoms. We explain the origin of this 
effect using elementary quantum mechanics, and summarize the experimental evidence for it. 
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I. INTRODUCTION 



Consider three identical bosons (or fermions with spin and isospin) that occupy a spatially 
symmetric s-state. Efimov predicted in 1970 that the spectrum obeys a geometrical scaling 
law, such that the ratio of the successive energy eigenvalues of the system is a constant. 1 
This scaling results in an accumulation of states near zero energy, with the size of the system 
growing larger by a factor of about 22.7 for identical bosons for every state approaching zero 
energy. The number of these three-body bound states is infinite when the dimer binding is 
zero. However, the number of three-body bound states is actually reduced as the two-body 
interaction is made more attractive. More remarkable is the fact that the Efimov effect is 
independent of the details of the interaction, for example, whether it is between atoms with 
a pair-wise van der Waals interaction, or between nuclei with a nuclear force. The effect 
persists even with unequal masses, as long as two of the three dimer bonds have near-zero 
binding. These statements will be made more precise in the following. 

A rigorous proof of Efimov's prediction followed in 1972.- Recently, an analytical solution 
for three identical bosons exhibiting the Efimov spectrum has been obtained.- Efimov exam- 
ined the triton and 12 C (the latter as a bound state of three alpha particles) for manifestations 
of his prediction. However, for nuclear systems, the interaction between two nucleons cannot 
be adjusted to yield a two-body zero-energy bound state. Also, the Coulomb potential be- 
tween protons often restricts the examples for which two of the three participants are loosely 
bound neutrons^ in the field of a nuclear core. However, for ultra-cold neutral atoms with 
tunable interaction strengths, these drawbacks are absent. The experimental observation of 
the Efimov effect was first made by Kraemer et aZ.- with an optically trapped dilute gas of 
133 Cs atoms at 10 nK. At such low temperatures, thermal motion does not mask quantum 
effects. Moreover, the two-body interaction between atoms may be fine-tuned using "Fes- 
hbach resonance. "— ~— Under such conditions, the recombination losses increase sharply due 
to the formation of Efimov trimers, giving the experimental signature of the effect. More 
recently, Barontini et al— have found evidence for two kinds of Efimov trimers in a mixture 
of 41 K and 87 Rb; namely KKRb and KRbRb. 

Although there are excellent reviews of Efimov physics (see, for example, Ref. lloL be- 
ginning graduate or senior undergraduate students unfamiliar with the quantum three-body 
problem might find them inaccessible. The Efimov effect arises from the large- distance 
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(asymptotic) behavior of the inverse square interaction. The spectrum resulting from such a 
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potential is well-known,— ~— and its implications for Efimov physics are discussed in Ref. 

The main goal of this article is to derive, using elementary quantum mechanics, the 
effective inverse square interaction in coordinate space of the three-body system. To achieve 
this goal, we use the separable form of two-body potential, for which the algebra is tractable. 
The Efimov effect is universal and does not depend on this special choice of a separable 
potential. A choice of the more familiar local form for the potential would entail more 
advanced aspects of the three-body problem, while not adding to the essential physics. Our 
hope is that this article can be used as the basis for introducing the fascinating Efimov effect 
in an advanced undergraduate, or first year graduate course in quantum mechanics. 

The outline of our paper is as follows. In Sec. II we briefly review the low-energy scat- 
tering of two particles.— The rationale for utilizing a nonlocal separable two-body potential 
near a resonance for the analysis of the three-body system^ is introduced at this stage. 
Short discussions of the separable potential, and Feshbach resonance follow. The relevant 
properties of the inverse square interaction for an understanding of the Efimov effect are 
given in Sec. III. Our discussion of the properties of the inverse square interaction sets the 
stage for the derivation of the Efimov effect in the three-body problem in Sec. IV. Recent 
experimental evidence for the Efimov effect, which has resulted in a renewed interest in this 
topic,— is briefly presented in Sec. V. Our concluding remarks are in Sec. VI. 



II. THE TWO-BODY PROBLEM 

A. Low energy parameters 

In the elastic scattering of two particles, the initial and the final energies are the same, 
and only momentum transfer may take place. The scattering amplitude is said to be "on the 
energy shell." The scattering cross section and other observables may be described in terms 
of a phase shift that occurs in each partial wave of the asymptotic (relative) two-body wave 
function. We assume that the two-body interaction is central, and falls off faster than r~ 2 , 
where r is the distance between the two particles. If three or more particles are involved, 
the situation is not so simple, even when the particles interact only pair-wise. If a third 
particle is present near the pair, it might take away some energy in the scattering process. 
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The two-body scattering amplitude is then said to be "off the energy shell," and depends 
on the short-distance behavior of the two-body wave function. 

The exact description of the three-body problem is given by three coupled integral equa- 
tions involving the two-body off-the-energy-shell scattering amplitudes.— The Efimov effect 
in the three-body problem is remarkable in the sense that its existence is dictated only by 
the low-energy two-body elastic scattering properties, and not by the off-shell behavior. The 
Efimov effect arises as long as the two-body interaction between at least two of the pairs 
is tuned to resonance, regardless of the interaction potential being, for example, local or 
nonlocal, of zero range or falling off as r -6 . At resonance, the two-body wave function is 
just about bound, and extends over a large distance. The tuning of the two-body interaction 
to resonance is all-important in giving rise to the Efimov effect. 

In this paper we choose a specific form of a nonlocal potential, called a separable potential, 
that gives an analytical form for the two-body wave function, and for which the treatment 
of the three-body Faddeev equations is simpler.- In the example of the three-body problem 
that we discuss in Sec. IV, we avoid solving the Faddeev equations explicitly, while still 
obtaining the Efimov spectrum. 

For two-body scattering at very low energies, the interaction is effective only in the 
relative s-state, because the centrifugal barrier in the higher partial waves suppresses the 
particles from coming close together. The s-wave scattering amplitude is well-known to be 
given by 



where 5o(k) is the s-wave phase shift. The differential cross-section for s-wave scatter- 
ing is isotropic and is given by \f (k)\ 2 . The scattering between the two particles is well- 
describecU^ 1 ^ at low energies by two shape-independent parameters, the scattering length 
a, and the effective range r , where 



The effective range r is a measure of the range of the two-body potential, and the 
scattering length is the intercept of the asymptotic zero-energy s-wave wave function. The 
reader unfamiliar with this description of the two-body problem should study Fig. 1 and 




(1) 



k cot S (k) = 1 — r k 2 + . . . 



(2) 



read Ref. QJ, Sec. II A. 
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Consider a situation where the attractive potential is not strong enough to support a 
two-body bound state, in which case the sign of the scattering length is negative.— As the 
potential is made more attractive, the scattering length becomes more and more negative, 
and goes to — oo at resonance (that is, a "bound state" at E — 0), and then flips to +00 
as the strength of the attractive potential is increased. A further increase in the attractive 
strength drives the scattering length to smaller positive values (but with a ^> ro), resulting 
in the formation of a dimer.™ For a zero-range potential, only the first term on the right- 
hand side of Eq. ([2]) survives, and it follows from Eq. ([I]) that the scattering amplitude is 
given by /o — — (l/a + ik)~ l . If we extrapolate the wave number to complex values k = ik,q 
for a bound state at k 2 = — Kq, we see that the resulting dimer has energy E b = -h 2 /Ma 2 , 
where M is the mass of each isolated atom. The corresponding wave function for r > r is 
given by 

, . . exp(— r/a) 

Mr) ~ 1 ' (3) 
r 

For r<a, Eq. ([3]) reduces to ipo ~ (1/r — 1/a), showing the connection of the intercept to 
the scattering length. 



B. Separable two-body potential 

A separable potential is a special form of nonlocal interaction which has been extensively 
used in the three-body problem in nuclear physics.— 1 ^ It results in considerable simplification 
of the Faddeev integral equation, which Amado and Noble 2 exploited to give a rigorous 
treatment of the Efimov effect. The work of Amado and Noble was further developed in 



Ref. 
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In the following we give the Schrodinger equation for two particles interacting by a 
separable interaction, and consider the scattering of a heavy particle of mass M interacting 
with a light one of mass m, with the mass ratio p = M/m. The relative energy is given 
by E2 = k 2 /u', where we have set h = 2m = 1 and v 1 = p/(p + 1). For simplicity, we 
assume p 3> 1, and take v' = \. The eigenvalue equation (after the center-of-mass has been 
removed) is 

(k 2 -H )\M=Vm, (4) 
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where Hq is the kinetic energy operator. In r-space Eq. fl4]) takes the form (where r is the 
relative coordinate between the two particles) 

(k 2 + V 2 )^ fc (r) = J (r|V|rV fc (r') dV, (5) 

and in p-space Eq. (jlj) becomes 

(k 2 - P 2 )Mp) = I (P\V\P')MP') d 3 p'. (6) 

A separable potential in Hilbert space may be written as V = — X\g)(g\, where the 
negative sign is taken for attraction, and A > determines the strength of the potential. 
In the coordinate representation, the attractive separable potential in the s-state is given 
by (r|V|r'} = —Xg(r)g(r'), where g(r) is taken to be real. When a bound state is present, 
the Fourier transform of g(r) is closely related to the bound state wave function. The 
Schrodinger equation in p-space for a bound state, k 2 = — Kq, is given by 



[4 + p 2 )^o(p) = Xg(p) I g(p')^ (p')d 3 p'- (7) 

In Eq. (J7D, 



g{p) = 

We may write Eq. ([7J) as 



J (p\r)g(r)d 3 r = ^ 3/2 J exp(-ip • r)g(r)d 3 r. (8) 



^(p) = AC ft0 .J^- v (9) 



g(p) 

(k 2 + p 2 

where C KQ = J ^(pO^oCp')^ 3 ^' ^ s a nonzero constant. If we multiply both sides of Eq. d^J) by 
g(p), and integrateover d 3 p, we obtain the equation that determines the binding energy Kq 
for a given potential 

g 2 {p) 



x !w^) dp=1 - (10) 

In the next section, a similar equation for the three-body bound state will be obtained. 
For explicit calculations we choose the popular Yamaguchi form 23 g(p) = (p 2 + /3 2 ) -1 , 

giving 

3 pKexp(-/3r) 



9( r ) = j^yjt J ex p( z P • r )g(p) d p = y 2 — r — • 

We substitute this choice of g(p) in Eq. (Q, and take its Fourier transform. After a little 
algebra, we find that 

, / x x ^ ft ( exp(-Kpr) exp(-/3r)\ 

Ao{r) = XC Ko y — ( I (12) 



For small binding, kq is only slightly greater than zero. We choose a short range potential, 
so that ^> kq, and asymptotically Eq. f fT2|) becomes the same as Eq. (J3J) obtained from 
the zero-range approximation. The universal nature of the Efimov effect is independent of 
the specific choice of g(r), as long as g(r) is short range. 



C. Feshbach resonance 



The strength of the interaction of a zero-range potential is controlled by the scattering 
length a. For cold-atoms the scattering length may be varied by a magnetic field, making 
use of the Feshbach resonance.- - - Near a Feshbach resonance, two unbound atoms with 
relative energy E slightly larger than zero may scatter off each other through a potential. 
This scattering process is said to occur in the "open channel." When the same two atoms 
come close together in a different total spin state, they may encounter a different potential, 
and form a bound or quasi-bound state at energy E res > E. The process leading to the 
formation of a bound or quasi-bound state is called a "closed" channel. The scattering 
state at energy E does not exist in the closed channel, because the closed channel requires 
E res > E to dissociate the atoms. Nevertheless, there may be interchannel coupling through 
a spin-dependent potential that might cause a virtual transition between the two channels. 
The two atoms in the open channel have different magnetic moments than when they are in 
the quasi-bound state, and hence the energy gap between the states may be controlled by the 
relative Zeeman shifts. As a result, the mixing between the two states is more pronounced 
when the Zeeman energy gap is reduced. In a single-channel description, 7 the two atoms in 
the open channel may make a virtual transition to the closed channel state and back, giving 
an effective single-channel scattering length in the open channel that is inversely proportional 
to the energy gap, AE, between the two states, where AE = Afi(B — B ). Here, A/i is 
the difference in magnetic moments of the two-atom system in the two channels, and B is 
the magnetic field at which the scattering length diverges (at resonance) and changes sign. 
When B is very different from the resonance value -Bo, the effect of the coupling is small, 
and the scattering length reduces to its natural "b ackg round" value abg. More details in the 



context of ultra-cold atoms may be found in Refs. 
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and 
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a < Energy a > 




FIG. 1: Adapted from Fig. 1 in Ref. LL6|. For a < 0, even though there is no two-body bound state 
(dimer), more and more three-body bound states (trimers) are formed as the potential becomes 
more attractive. The solid lines represent these Efimov trimer states which emerge at E = 0, and 
a < 0. For a > on the right, the trimers break up at the atom-dimer continuum , indicated by 
the dotted curve 

III. EFIMOV SPECTRUM AND THE INVERSE SQUARE POTENTIAL 

For three identical bosons interacting pairwise via short-range potentials, Efimov pre- 
dicted an infinite number of three-body bound states with geometrical scaling, when the 
two-body interactions are resonant, that is, \a\ — > oo. 



Figure 1 (following Ref. \lw depicts the Efimov scenario where the energies of some of the 
three-boson states are plotted as a function of 1/a; for E > the energies of the three-atoms 
form a continuum. The bound Efimov states are shown schematically by solid lines, with a 
scaling factor (ratio of successive energy eigenvalues) set artificially at 2 rather than 22.7. 
The Efimov states break up on the positive side of a in the atom-dimer continuum given by 
-h 2 /(Ma 2 ) < E < (see the dotted curve in Fig. 1). 

The Efimov spectrum of the three-body problem is the signature of an attractive central 
potential which falls off asymptotically as the inverse square power of the distance. A three- 
body long-range potential is suggested by noting that for large positive a, the size of the 
dimer is very large, and the presence of another atom, even if very distant, may be "sensed" 
by the dimer. For large negative a the two atoms, even if not bound, are spatially correlated 
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over a distance of order \a\ in a quasi-bound state. 

To obtain the geometric scaling of the spectrum, it is sufficient to consider the simpler 
problem of a single particle of mass m in an inverse square potential V(r) = (h 2 /2m)\/r 2 , 
where A is a dimensionless coupling constant. Classically the equation of motion in this 
potential is scale invariant under the continuous transformations r — » ar, and t — > a 2 t. 
Quantum mechanically, for A > —1/4, there is no bound state, and the continuous scale- 
invariance is valid. A zero-energy state appears for A = —1/4, and the system is anomalous 
for A < — 1/4, due to the short- distance singularity of the potential.— 1 ^ A direct consequence 
of the anomaly for A < —1/4 is that there is no longer a lower limit in the energy spectrum, 
and a regularization is required.— We are interested in the situation where the potential is 
inverse square only for r > r c , where r c is taken as the short- distance cut-off. We impose 
the boundary condition that the eigenfunctions vanish at r = r c , which results in a discrete 
spectrum. The geometric scaling property, namely that the ratios of the adjacent energy 
eigenvalues remain a constant, is independent of r c . 

We write the Schrodinger equation in the s-state for r > r c as 



d 2 (s 2 + 1/4) 



2 772 

u(r) = -^Eu(r), (13) 



dr 2 r 2 

where at this stage, Sq > is just a way of parametrizing the strength of an inverse square 
potential that is greater than 1/4. For bound states, we set (2m/h 2 )E = — k 2 , and require 
that wave functions vanish at infinity. We then obtain the solution u(nr) = y/KrKi So {nr), 
where Ki SQ is the modified Bessel function of the third kind of pure imaginary order iso.^ 
The boundary condition that u(nr c ) = makes k discrete, such that K iSo (K n r c ) = 0, with 
n a positive integer. For shallow bound states such that (ft„r c 1), the zeros of the Bessel 
function K iso (K n r c ) are given by 

n n r c = exp (— )(2 e ^)[l + O(s ) + ...], (14) 

V So / 

where 7 is Euler's constant. Equation (14) leads to the desired result 



exp(— 2n/s ), n — 1,2, ...00 (15) 



which is the geometric scaling mentioned previously. Note that the actual value of E n scales 
as r~ 2 , but geometric scaling holds for the shallow states. Also note that as n becomes 
larger, the states become shallower, with an infinite number of states accumulating at zero 



energy. In the three-body problem in which three particles interact pair-wise, there are six 
degrees of freedom after the center-of-mass motion is eliminated. This problem is commonly 
treated in hyperspherical coordinates^ with a hyperradial variable R, and five angles. For 
equal mass particles,— R = \J (rf 2 + r% 3 + r| 1 )/3. In the adiabatic approximation^! for fixed 
R, we solve the Schrodinger equation for the angular variables, thus obtaining a complete set 
of adiabatic eigenvalues s(R) and corresponding eigenfunctions. The solution shows that in 
the resonant limit a — > ±00, neglecting channel coupling, the same inverse square potential 
as in Eq. (1) appears in hyperspherical coordinates, with r replaced by R. 

For identical bosons, exp(7r/s ) — 22.694, (that is, Sq ~ 1.00624) but in general, Sq 
depends on the mass ratios. If a is finite and very large, with \a\ 3> r , the inverse square 
interaction e(R) is cut-off at a short distance of the order of tq, and at a long distance of 
the order of \a\. The number of shallow bound states is given approximately by 15 

iV-^hiH (16) 
7T r 

A simplified derivation of Eq. f fl6|) is given in the Appendix. 



IV. THE THREE-BODY MODEL 



Our objective of obtaining an inverse square interaction in the three-body problem is best 
served by taking two identical heavy particles 1 and 2, each of mass M, and particle 3 of 
mass m <C M, as shown in Fig. 2. 

The coordinates of the particles are labeled r 1; r 2 , and r 3 , as measured from an arbitrary 
origin not shown in Fig. 2. The three-body analysis is simplified if the relative coordinates 
R = (ri — r2), and r = r3 — (ri + r 2 )/2 are introduced. Following the common convention, 
we denote by V\ the interparticle potential between particles 2 and 3, and likewise for V2 
and V3 in cyclic order. The three-body Schrodinger equation is given by (p = M/m, and 
recall that h = 2m = 1) 

fftf(r,R) = £tf(r,R), (17) 

with 

H = ~V 2 R -^V 2 r + V 1 + V 2 + V 3} /i = p/2, u = 2p/(2p+l), (18) 

where E is the energy of the three-body system. For p ^> 1, v — > 1, and the motion of 
the heavy particles is very slow compared to that of the light particle of mass m, which is 
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FIG. 2: A schematic illustration of the relative coordinates used in the three-body problem. See 
the text for details. 



in the spirit of the Born-Oppenheimer approximation. Following Ref. [15| we may apply the 
adiabatic Born-Oppenheimer approximation to solve the three-body Schoodinger equation 
in two stages. First, the wave function is decomposed according to 

*(r,R) = V(r,RMR), (19) 

where ?/>(r,R), with eigenenergy e(R), is first solved for the relative motion of the light- 
heavy system, keeping R a parameter. For fixed R, the relative kinetic energy of the heavy 
particles is zero, and the potential energy V 3 (R) is a constant shift in energy. Equation ffTTj) . 
with v — 1, becomes 

[-V* + V 1 (r - R/2) + V 2 (r + R/2)] ^(r, R) = e(R)^(r, R). (20) 

Next, the two-body Schrodinger equation with e(R) as the adiabatic potential between the 
two heavy particles is solved for the ground state energy, E, within the adiabatic approxi- 
mation: 

[-V|//i + V 3 (R) + e(R)] <f>(R) = E<P(R). (21) 

Because the interatomic potential V 3 (R) falls off faster than 1/R 2 , it does not affect the 
asymptotic behavior of e(R), and we set V 3 (R) = 0. For V\ and V 2 in Eq. fl20|) . we take 
short-range separable potentials — X\g)(g\. Thus, for a fixed parameter R, we have 

(r' - R/2|V!|r - R/2) = -A(r' - K/2\g)(g\r - R/2), (22) 

11 



and likewise for V2, with — R/2 replaced by +R/2. 

If we introduce the displacement operator, D R / 2 = exp (ip r - R/2), we can write |r + 

R/2) = Dji/ 2 \t), where p r is the momentum conjugate to r. In operator form, Eq. (l2"0|) 
becomes 

[-V? - A (D\g)(g\D- 1 + D~ l \g){g\D)~\ |^> = e (J2)|V>), (23) 

with the subscripts on D suppressed for brevity. Equation ( 1231) bears a striking similarity 
to the two-body equation in operator notation given by [— V 2 . — A|g)(g|] \ip) = E 2 \ip). In 
Eq. (J23D, let e(R) = -k 2 (R), G = (-V 2 + k 2 )~\ and (g\D\ip) = N(R/2). The light-heavy 
wave function ?/>(r, R) has a definite parity under core-exchange, R — > — R, so that for the 
lowest energy state with energy E [recall that we are solving Eq. (1211) ]. we take the wave 
function to be symmetric. We define (glD^ 1 ^) = iV(— R/2) = iV(R/2), thereby reducing 
Eq. ([23D to 

\^} = XG(D- 1 \g} + D\g})N(R/2). (24) 
We multiply both sides by (g\D, note that [G, D] = and iV(R/2) 7^ 0, and obtain 



y f (g\g) + (g\D 2 \g) 
V (-v r 2 + K 2 ) 



1. (25) 



In the momentum representation, Eq. (20) takes the form 

9 2 (Pr) ,3 , \ f 9 2 (Pr) exp(ip r ■ R 



J pf + K 2 J p 2 + K 2 

which is analogous to Eq. ( ITU]) for the two-body problem. Note that Eq. (l2"oT) assumes p 1, 
so that v' ~ z/. Recall that v' = p/(p + 1) and z/ = 2p/(2p + 1). The coupling constant A 
may be eliminated by fixing the binding of the two-body problem. For the Yamaguchi 
form g(p) = (p 2 + /3 2 ) -1 , the integrals in Eq. (|26|) may be performed analytically, giving 



(27) 



(/3 -k) 2 R /3 -k 

As the distance R between the two heavy atoms is increased, the light atom will tend 
to attach to one of them, and k 2 (R) — > Kq. We are particularly interested in this large R 
behavior of e(R) as the two-body binding Kq goes to zero. To deduce the functional form 
of e(R) = —k 2 (R) for large R, it is convenient to define £ = k — Kq, and substitute it into 
Eq. ( 1271) . From Eq. (llip . we see that g(r) is short-ranged for /3 > 1. In the resonant limit, 
consider a — > 00 from the positive side. Then k = 1/a — > 0. We consider two possible cases 
in turn: 
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(1) k < 0, 0R > 1, R/a -> 0, £ < /3 for which exp(-/3i?) -» 0, and exp(-K -R) -> 1- 
Some algebra yields the equation 

P -£R 

W " l ' (28) 

which has the solution £R = A, where A = 0.5671 . . .. Hence k = kq + £ = kq + A/R. 
In the limit of a — > oo, Ko — > 0, we see that e(i2) = — k 2 = —A 2 /R 2 , giving the desired 
inverse square potential. In our analysis, we set v = 1, which only approximately holds 
for p 3> 1. However, setting v = 1 is not necessary, and Fonseca e£ a/.— have shown that 
e(R) = —A 2 /{yR 2 ). The effective potential including the kinematic factors that appear in 
the equation analogous to Eq. (fT3"j) may be easily deduced from Eq. f[2~Tj) to be — pA 2 / [vR 2 ) = 
— 1/AA 2 (1 + 2p)/R 2 . Therefore the scaling parameter so in the Efimov spectrum, as defined 
in Eq. (Tl3|) . depends on the mass ratio p. 

(2) Let the scattering length be large but finite, and consider a distance R that is even 
larger, so that R/a ^> 1. R/a ^> 1 implies that kqR 3> 1, even though the other conditions 
are the same as in case (1). Thus, in spite of being able to set exp(— (3R) = 0, we cannot 
set exp(— KqR) to unity. We may show that for R 3> a, Eq. (12"7|) yields the Yukawa form 

o p-R/a 

This form of a static Yukawa potential arises naturally in nuclear physics due to the exchange 
of a light mass boson (for example, a pion) between two heavy mass nucleons. In nuclear 
physics the range of the potential is determined by the square root of the mass of the light 
particle; in our case, it is the square root of the binding energy of the light atom that plays 
the analogous role. This situation is depicted in Fig. 3. In region I, R < Rq, where Rq 
is the range of the interatomic potential 2 ^ between the two heavy particles. We cutoff the 
potential for the shallow states at this distance. In region II there is a 1/R 2 potential, 
which extends to distances R a, hence for all R as a — > oo, as in case (1). Region III 
is for R ~ a, where the transition to the Yukawa form takes place. Region IV contains 
the asymptotic behavior of the potential. As the two-body binding Kq goes to zero, the 
long-range potential e(R) takes the inverse square form, which has no length scale. From 
dimensional considerations, because h 2 /M has the dimensions of EL 2 , an inverse square 
potential is the only form possible in the absence of other mass scales or coupling constants. 

From Fig. 3, a related prediction associated with the collapse of the three-body system 
called the Thomas effect^ may also be deduced. Consider a two-body system of range r 
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FIG. 3: The effective potential s(R) between the two heavy- mass particles that arises in the 
adiabatic approximation due to the interaction with the light-mass particle. For R < Rq, the 
interatomic potential is irrelevant for the Efimov effect. The scattering length a is marked on the 
horizontal axis. As a — > oo, the two-body binding Kq goes to zero, and e(R) — > R~ 2 for all R > Rq. 
The four regions of the potential e(R), are discussed in the text. 

with a fixed binding Kq. Let the range parameter r become smaller and smaller, adjusting 
the strength of the two-body potential so that k$ remains constant. The Thomas effect 
asserts that the three-body system will collapse as tq — > 0, with its deepest bound state 
going to — oo. In Fig. 3, the short distance cut-off i? of the inverse square potential e(R) 
goes to zero as r goes to zero. Such a behavior near the origin causes a collapse in the three- 
body energy. The Thomas effect does not require the asymptotic form of the potential to be 
inverse square, but rather is associated with the short distance singularity of the potential. 
From Eq. (fTrJj) we note that the number of three-body bound states diverges when the ratio 
a / r o oo, which may be brought about either by letting a — > oo with ro finite, as in the 
Efimov effect, or letting r — > 0, with a finite, which corresponds to the Thomas effect .~ 
Unlike the Efimov effect, the Thomas effect is not amenable to experimental verification 
because the range of the two-body potential cannot, as of now, be tuned to zero. 
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V. EXPERIMENTAL EVIDENCE 

Starting with the pioneering work of Kraemer et a/.- with ultra-cold Cs atoms in 2006, 
several experiments^— have confirmed Efimov's predictions by measuring the three-body 
recombination losses of atoms through the reaction A + A + A — > A 2 + A. The experiment 
with heteronuclear atoms was done by Barontini et al 9 All these experiments require very 
low temperature fcgT < h 2 /Ma 2 for a very large scattering length a to avoid break up of the 
dimers. For example, for Cs atoms, T = 10 nK was required to largely eliminate thermal 
effects. 

The atoms in a gas have a tendency to form lower energy dimers directly for scattering 
length a > 0. But two atoms alone cannot form a bound state and preserve momentum 
and energy at the same time. Dimer formation is possible if a third atom is nearby within a 
distance of the order a. This problem was studied in a gas of identical atoms with number 
density n, while looking at atomic losses in Bose-Einstein condensates.-^ Let the number 
of three-body recombinations per unit volume per unit time be denoted by i/ rec , which is 
proportional to n 2 (av)(na 3 ). Here n 2 (av) is the probability of two atoms being in the 
interaction volume (erf), where a is the elastic scattering cross section and v is the relative 
speed between the two atoms. The probability of finding a third atom within a distance a 
is na 3 . Because v = hk/m ~ h/(ma), and o oc a 2 , we obtain 



where m is the mass of each atom, and C(a) is a dimensionless constant. The variation 
of C(a) with the scattering length a exhibits the emergence of an Efimov trimer (negative 
a), and its subsequent dissolution to dimer plus atom (positive a), see Fig. 1. It has been 



mass of the trimer, and the fact that three atoms are lost for each trimer. To understand the 
experimental peak in the measured recombination length p3 (see Fig. 4), note that shallow 
dimers are not formed directly for a < 0. The potential e(R) of Fig. 3 develops a potential 
barrier for large R, followed by an attractive potential at shorter distances. Such a potential 
emerges from a hyperspherical calculation for equal mass bosons.— There may result, for 
certain negative values of a, shape resonances of the trimer at E > 0. When the energy of the 
three atoms matches a resonance, there is enhanced barrier penetrability to form an Efimov 




(30) 



calculated from theory.— — The recombination length (the ordinate of Fig.4) is defined as 
P3 = (2v / 3C(a)) 1//4 a, where the numerical factor is included to take into account the reduced 
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FIG. 4: Adapted from Fig. 2 of Ref. I37l . The solid circles are experimental data points^ at 10 nK 
for the observed three-body recombination: Cs + Cs + Cs — > C S2 + C s. The scattering length is 
measured in units of the Bohr radius ao • 

trimer. The enhanced barrier penetrability allows the system to relax to a deeply bound 
dimer state (from the interatomic potential), with the excess energy being carried away by 
the third atom. The peak on the negative side is at a = — 850ao, where ao is the Bohr 
radius. The next peak should be at a value of a that is 22.7 times larger, that is, a = —19, 295, 
which is outside the range of experimental observation. Nevertheless, these multiple peaks 
were observed in later experiments.— ~— Even though multiple resonance peaks were not seen 
in the Kraemer et a/.- experiment, the interference minimum observed on the positive a side 
is convincing evidence for the Efimov effect.— Hyperspherical calculations for a > not only 
give an attractive potential as in Fig. 3 for s(R), but also a barrier for positive energies that 
falls off for large R. The three incoming atoms may be reflected by the barrier, and form 
A 2 + A on the way out, interfering with the incoming dimer plus atom system following the 
large- R attractive channel of Fig. 3. The out-of-phase interference between the two paths 



leads to the curve that fits the experiment for positive a (see Fig. 2 of Ref. 1371 ) 
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VI. CLOSING REMARKS 



It is remarkable and not foreseeable in 1970 that the Efimov effect has been experimentally 
verified using ultra-cold atoms. We have considered a simple case for which two of the atoms 
are very heavy compared to the third. The light heavy pairs are treated very differently, as 
evidenced by the adiabatic approximation. The light particle is almost unbound, and can 
be very distant from the heavy partners. The overall size of the shallow Efimov trimer is 
large in comparison to the short-range of the two-body interaction. 

As is clear from Sec. IV, it is the dynamics of the light particle that generates the at- 
tractive long-range interaction between the heavy ones. We may regard the light particle as 
being exchanged back and forth to generate this potential between the two heavy particles. 
Even for equal mass particles, Efimov^ has stressed this intuitive picture for the inverse 
square potential. A different approach must be taken for three identical particles, as indi- 
cated briefly in Sec. Ill for three identical bosons. Each particle has to be treated the same 
way, and hyperspherical coordinates are well-suited to describe the shallow Efimov states. 
The size of these states is large, and the small binding results in the formation of floppy 
triangles. For these shallow states, the coupling between the hyperradial variable R and the 
hyperangles is small. The coupling vanishes for \a\ — > oo, and the hyperangular contribu- 
tion just adds a centrifugal term for angular motion to the adiabatic potential V(R). The 
deepest potential is attractive enough (for L = 0) to yield the Efimov effect. 

For three nucleons the overall wave function under the exchange of any two has to be 
antisymmetric. As was originally found by Efimov, 1 imposing these symmetries results 
in only the spatially symmetric angular momentum state L = (with the spin-isospin 
combination antisymmetric) leading to an attractive interaction. Unlike charge neutral ultra- 
cold atoms, it is not possible to manipulate the strength of the nucleon-nucleon interaction. 

Given that there is suggestive experimental evidence for the Efimov effect in the four- 
body sector (see for example, Ref. 133), it is clear that Efimov physics in ultra-cold atoms 
will continue to be an active area of research. 
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Appendix: Derivation of bound states in the inverse square potential 



One way to obtain the number of bound states of a potential is to calculate the canonical 
partition function Z(/3) = ^ exp(—(3Ei), where (3 for our purpose may be taken to be a pos- 
itive parameter. The expression for Z(f3) may be written as Z{j5) = J °° g(E) exp(—(3E)dE, 
where g{E) = 2~2i$(E — Ej) is the density of states. The latter equation shows that Z((3) 
is the Laplace transform of g(E). The inverse Laplace transform of Z((3) with respect to E 
yields the density of states. The integration of the density of states in an energy interval 
gives the number of states. In the following, we present a simple semiclassical derivation of 
the number of shallow Efimov states in the energy interval between and E. 

Consider the one-body inverse square potential given in Eq. (1) 

;2 /„2 



v{r) = _pA±m. (A.l) 
2m r z 



This form of the potential is taken to be valid for r$ < r < a and for large a (it is assumed 
that ro — r). We shall now derive Eq. (3). The semiclassical partition function for a given 
partial wave is 



1 



Ze(P) = il dp r exp(-l3pi/2m) / drexp(-^(r)). (A.2) 



z t (p) = ( ) /3-i/2 / e -m(r) dr _ (A 3) 



After performing the momentum integral in Eq. f ]A.2j) . we obtain 

\ 1/2 ^oo 

,2vr^ 

The Laplace inversion of Zg(/3) with respect to f3 gives the density of states 

s ' {E) -{w)\S^Sm e(E - v ' [r)y (kA) 

For the inverse square interaction a WKB-type approximation yields exact results when the 
Langer correction is implemented,— that is, £(£ + 1) is replaced by (£ + 1/2) 2 . Hence, 



ft 2 (4 + 1/4) , 

2m r 2 2mr 2 



Vi(r) = -— v ' J + ^—A£ + 1/2Y. (A.5) 



We are interested in the I = partial wave, for which 
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The number of states is obtained by integrating Eq. (1A.4j) with respect to E, followed by 
the integration with respect to r: 

NiE ^{w)l£^^ dr - (A7) 

If we take E close to zero and Vo(r) from Eq. (\A.6\ we obtain 



^ln(^), (A.8) 
7i \r n / 



N 

>o 

which is the well known result for the number of shallow Efimov states.- The same derivation 
applies to the three-body problem with potential e(R) (see Fig. 3), when, for very large a, 
we take the form Eq. flA.ip with Rq < R < a. 
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